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We show that double-layer graphene (DLG), where an external potential induces a charge- 
imbalance between n- and p-type layers, is a promising candidate to realize an exciton condensate 
in equilibrium. To prove this phenomenon experimentally, we suggest coupling two DLG systems, 
separated by a thin insulating barrier, and measuring the excitonic Josephson effect. For this pur¬ 
pose we calculate the ac and dc Josephson currents induced by tunneling excitons and show that 
the former only occurs when the gate potentials of the DLG systems differ, irrespective of the phase 
relationship of their excitonic order parameters. A dc Josephson current develops if a finite order- 
parameter phase difference exists between two coupled DLG systems with identical gate potentials. 


The search for the long ago predicted excitonic insu¬ 
lator (El) state has recently stimulated a lot of experi¬ 
mental work, e.g., on pressure sensitive rare-earth chal- 
gogenides, transition metal dichalcogenides or tantalum 
chalcogenides^^®. Theoretically the excitonic instability 
is expected to happen, when semimetals with very small 
band overlap or semiconductors with very small band gap 
are cooled to very low temperatures®’^. To date there ex¬ 
ists no free of doubt realization of the El, however, and 
even the applicability of the original El scenario to the 
above material classes is a controversial issue®’®^^®. There 
are serious arguments why the El in these bulk materi¬ 
als, if present at all, resembles rather a charge-density- 
wave state than a “true” superfluid exciton condensate 
exhibiting off-diagonal long-range order^^’^^. 

On these grounds a non-ambiguous experimental proof 
of a macroscopic phase coherent exciton condensate 
would be highly desirable. Spectroscopic analyses have 
not established an exciton condensate so far. The char¬ 
acteristics of junction devices, where at least in one com¬ 
ponent an El is realized, may lead to valuable insights in 
this respect^®. Due to the proximity effect a high resis¬ 
tance should appear across a semimetal-El junction that 
distinctly differs from that of a semimetal-semiconductor 
device^^. In coupled quantum wells, Josephson oscilla¬ 
tions should accompany exciton condensation^®’^®. Here 
we will pursue a similar idea, namely that a Josephson- 
type tunnel current might appear when two El systems 
are coupled to each other by a thin insulating barrier such 
that coherence is established between the condensates. 

Two-layer systems of spatially separated electrons and 
holes that feature an attractive interlayer electron-hole 
coupling are particularly suitable for a Josephson-type 
tunnel experiment. In this case a condensate of excitons 
might occur when the tunneling between the layers is 
negligible, but the corresponding Coulomb interaction is 
not^^. Double-layer systems thereby inhibit the obsta¬ 
cles coming from interband transitions or the coupling to 
phonons, which inevitably occur in bulk materials and 
prevent a possible exciton condensation by destroying the 
U{1) symmetry^^’^®’^®. It is also advantageous that (ex¬ 
citon) tunneling effects are experimentally well accessible 


in double-layer systems. 


Ix^ 





FIG. 1. (Golor online) Upper panel: Schematic setup to ver¬ 
ify an excitonic Josephson effect. Two DLG systems are sepa¬ 
rated by a thin tunnel barrier (TB). By external potentials Vi 
and Vr the chemical potential (particle numbers) can be tuned 
in each layer of the left and right DLG, respectively. Here the 
charge carriers are electrons (e) and holes (h) in the upper 
(lower) n-type (p-type) layer. The layers are separated by a 
dielectric of thickness d. Tunneling of coherent electron-hole 
pairs (excitons) will induce an electron current Ix through 
the barrier. Obviously, the current in the lower layer equals 
the current in the upper layers in modulus but Hows in the 
opposite direction. Lower panel: Band structure near the K 
point and chemical potential fii of neutral DLG with Vi — 0 
(left) and gated DLG where 14 > 0 (middle). First Bril- 
louin zone of DLG (red hexagon) with high-symmetry points 
(right). Backfolding parts I and II, the grey rectangle shows 
the representation of the Brillouin zone we use in Fig. 2. 

Graphene-based double layers (separated by an ade¬ 
quate dielectric, e.g., hexagonal boron nitride or Si 02 ) 
show great promise for realizing a corresponding setup 
(see Fig. 1). For double-layer graphene (DLG), a gate- 
bias across the layers creates a charge imbalance, where¬ 
upon the attractive Coulomb interaction between the ex¬ 
cess electrons and holes on opposite layers raises the pos- 
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sibility of exciton formation. While a fine tuning of the 
band gap can be achieved by the external potential, the 
recombination of electrons and holes can be fully sup¬ 
pressed by the dielectric^®. Then, in the weak coupling 
regime, exciton condensation is triggered by a Cooper- 
type instability, where the particle-hole symmetry of the 
system ensures a perfect nesting between the electron 
Fermi surface and its hole counterpart in the n- and p- 
type layers of a biased DLG system^^^^®. Placing an (in¬ 
sulating) tunnel barrier between two such DLG systems, 
the Josephson current can be used to analyze whether or 
not an exciton condensate has been formed in each of the 
subsystems. We note that a similar setup was recently 
proposed to study thermal transport in a temperature- 
biased exciton-condensate junction^"*. 

The tight-binding Hamiltonian we assume for a DLG 
subsystem i (i = l,r, left or right) has the form 

k k 

+ ]^ ^ ^k,k',q Ok+qiOki^k'-qi^k'i > (1) 

k,k',q 


where and annihilate (create) electron quasipar¬ 
ticles in the n layer and p layer, respectively, with in¬ 
plane momenta k and band dispersions 


= ±70 


3 -I- 2 cos (^VSky^ 

-1-4 cos {Viky/‘2}j cos {3kx/2) 


T ■ (2) 


In the low carrier density regime, the effective band struc¬ 
ture ( 2 ) should account for the effects of the intralayer 
Coulomb interaction. The corresponding particle trans¬ 
fer amplitude is parametrized by 70 — 2.8 eV (which 
defines the unit of energy in what follows) and the mo¬ 
menta kx,y by a~^, where a ~ 0.142 nm is the carbon- 
carbon distance within graphene’s honeycomb structure 
(see Ref. 25). The external potential Vi determines the 
chemical potential: pi = Vi/2. The interlayer Coulomb 
interaction leading to exciton formation is 


C^k,k',q 


g-<i|q| 

= K —;— COS 

|q| 




( 3 ) 


where k, = gs2TT/e, = 2, e denotes the dielectric con¬ 
stant of the dielectric, and N gives the total number of 
particles^®’^^. Since electron-hole recombination is pre¬ 
vented by the dielectric, we neglected in Eq. (1) all in¬ 
terlayer Coulomb interaction terms that do not preserve 
the number of electrons (or holes) in a single layer, e.g., 
Hu oc aL+q.akiaL-q^V* or Hu oc al,+q,&kiaL-q»V*- 
Note that the model (1) exhibits a U{1) symmetry, which 
causes the phase of the El order parameter to be unde¬ 
termined. 

A mean-field decoupling of the Coulomb interaction 


yields 


Hi — ^ ^ki^L^ki + ^ ^ki^L^ki 

k k 

± ^ ^ki^L®ki + ^ki^kj^ki’ (4) 

k k 

where we have introduced the El order-parameter func¬ 
tion 
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where = ©k+q - ©k (4'' = ©k'+q - ©k') with ©k = 
arctan(A:y/fca;). The phase of A^- determines the phase 
of the ground-state wave function. 

In view of the specifics of the graphene spectrum, 
which are suggestive of large screening effects, the cor¬ 
rect approximation for the screening of the interlayer 
Coulomb interaction has been a controversial issue^®^®®. 
A mean-field inclusion in the normal phase certainly over¬ 
estimates screening and reduces |Ak| in an unrealistic 
way®®. In addition, correlations between electrons and 
holes may substantially weaken the screening®®. As a 
result double-layer graphene-based systems are not that 
disadvantageous for the realization of an excitonic con¬ 
densate as naively might be expected. Determining the 
transition temperature of the El phase necessitates, of 
course, an appropriate treatment of (dynamical) screen¬ 
ing in the condensed phase^®’®®’®®. This is beyond the 
scope of this work. In order to discuss the excitonic 
Josephson effect in the case of a realized El ground state 
at zero temperature, Eqs. (I)-(5) are adequate to leading 
order. 

We first solve the self-consistency equations for the El 
order parameter at zero temperature, assuming k = 7.0 
and d = 2.5. Figure 2 shows its modulus |Aki| within the 
first Brillouin zone. Apparently only electrons and holes 
near the Brillouin zone’s K or K' points are bound into 
excitons. These particles occupy the states closest to the 
Fermi energy and therefore will be most susceptible for 
electron-hole pairing. Since the external potential fixes 
the position of the Fermi energy, it determines the behav¬ 
ior of Aki as well. When V) is raised from zero, more and 
more states become available for an electron-hole pairing. 
As a result the El order parameter function is finite in a 
larger region of the Brillouin zone and the total (momen¬ 
tum accumulated) order parameter increases. At Id « 2 
the El order parameter attains its maximal value and 
starts to decrease if the gate voltage gets larger until the 
El phase breaks down at about 1^ « 6 , where the Fermi 
energy coincides with the upper (lower) edge of the con¬ 
duction (valence) band. The relatively sharp boundary 
confining the area of bound electron-hole pairs indicates 
a BCS-like pairing®^. 

A Josephson effect occurs when two DLG systems, 
where quantum coherence is realized, will be coupled 
to each other®®. Assuming that the ground states in 
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FIG. 2. (Color online) Modulus of the El order parameter 
function for DLG in reciprocal space, |Aki|, at different gate 
voltages Vi. The upper line plots give |Aki| along the high- 
symmetry directions of the Brillouin zone (cf. Fig. 1), for (a) 
Vi = 0.0 (black dotted line), Vi = 0.5 (blue dashed line) and 
(b) Vi = 1.0 (green dot-dashed line), Vi = 2.0 (red solid line). 
Corresponding intensity plots in the first Brillouin zone (using 
a 20 X 20 grid in k-space), for (c) Vi = 0.0 and (d) Vi = 2.0. 


both subsystems are described by macroscopic wave func¬ 
tions, {i = l,r index the left, respectively, 

right, subsystem), which—for simplicity reasons—are as¬ 
sumed to be equal in modulus but may have a different 
phases. The current density that flows between the left 
and right subsystem is j = ^ (T'jV'F* — (we set 

both h = 1 and the electron mass Wg = 1; e denotes the 
electron charge). Obviously any finite phase difference 
A(j) = (j)r — 4>i prevents the current density from vanish¬ 
ing, i.e., a persistent tunnel current flows through the 
barrier. 

The Hamiltonian of the coupled system is H = Hi + 
Hr + Ht , with 

Ht — Tk,p (^ahupr + 

k,p 

+ ^k.p (aprOki + ^pr^k/) > (6) 

k.p 

describing the tunnel process. 

Tunneling excitons cause an electron current in the 
n layer. This current equals the one in the p layer— 
which flows in the opposite direction, however—in mod¬ 
ulus. Below we adapt the approach outlined in Ref. 34 
to the (coherent) exciton tunnel processes in DLG. For 
this we define the tunnel current as the time derivative 
of the number of electrons in the upper layer and use an 
^-matrix expansion approach^^. 


I(t) = e{Nai{t)) = -ie 


d£ 


Nal{t),HT{t') ) , (7) 


where Nai = Ek4iaki> = -* Ek.p ^k.paL^pr + 
*Ek.p?i:.p4^ako and 


Hrit') = , with H' = Hi + Hr. Note 

that the chemical potential in the left and right DLG 
systems may differ. We introduce the applied junction 
voltage 


W = pr- Pi, ( 8 ) 

and the operators A{t) = X^k p ^k,p®pr(^)“ki(^) 

B{t) = Ek.p^k,p^pr(0^ki(^)- Although the specific 
choice of the tunnel matrix element will affect the Joseph- 
son current, we leave the investigation of this quantita¬ 
tive facet for future work. Instead we focus on the influ¬ 
ence of different chemical potentials and phases of the El 
order parameters in the DLG systems, which will domi¬ 
nate the physics. Analyzing Josephson tunneling for nor¬ 
mal superconductors the combination Tk^pT^ p is usually 
approximated by |Tp times a phase factor, where both 
quantities are assumed to be independent of k and 
In our case, the order parameter function has a rather 
complex momentum dependence, however (see Fig. 2). 
To avoid a fourfold numerical integration over the Bril¬ 
louin zone, we employ Tk.p = ^k.p, which—in the sense 
of a “principle of proof” calculation—will underestimate 
the effect we are looking for. Any more sophisticated 
treatment of the tunnel barrier will increase both the 
tunnel and normal currents. 

With that, the tunnel current takes the form 

m = -{ife f dt' (([A(t),At(t')])e*'^(*-‘') 

+ {[A{t), B^t')]) + {A{t), A{t')]) 

+ ([A(t),H(t')])e*^(‘-*') - ([At(t),A^(f')])e-*’^(*+*') 
- {[A^{t), - {[A\t), A{t')]) 

-([At(t),H(t')])e-*'^(‘+‘')). (9) 

Let us first focus on the current due to tunneling exci¬ 
tons exclusively. This current is governed by the second 
term and eighth term of Eq. (9): 

Ixit) = 2elm [e^^^^Xreti-W)] , (10) 

containing the retarded Green’s function Xret{W) = 

- O, Xretit-f) = 

t'){[A{t),B^{t')]). 

Normal tunneling electrons are described by the first 
and the seventh term of Eq. (9): 

/„(t) = -2eImrret(VF), (II) 

where Wet(fE) = = 

—i&{t — f) ([A(t),A'l'(t')]). It is worth noting that split¬ 
ting up the current into a quasiparticle (normal) and in¬ 
terference (oscillating Josephson) current was also carried 
out analyzing the thermal transport through Josephson 
junctions based on DLG^"*. 

We now factorize X^etiW) and Tret(kF) into contribu¬ 
tions stemming from the left and right subsystems. For 










4 


these we use the mean-field Green’s functions calculated 
with the Hamiltonian (4), which are expressed in terms 
of the quasiparticle states. 

Figure 3 gives the time dependence of the tunnel cur¬ 
rent for four characteristic situations. We first consider 
the case that two identical DLG systems are coupled by 
a thin barrier, i.e., Vi = Vr, W = 0. Clearly, if both El 
order parameters have the same phase, no current flows 
through the junction. A dc current arises when the left 
and the right systems differ in terms of the phase of the El 
order parameter: A(j) = (j)^ — (pi = —7r/2 in Fig. 3. An ac 
current appears, on the other hand, if a finite (constant) 
voltage is applied across the junction. Its frequency is 
2W, which coincides with the frequency of the Joseph- 
son current for coupled superconductors. An additional 
phase difference amplifies the tunnel-current amplitude 
and leads to a phase shift in the current. 

We finally analyze the voltage-current characteristics 
as to its phase dependence. For this purpose, we switch 
off the gate voltage in the left DLG system (V/ = 0) and 
fix its El phase to tt. Now the gate voltage in the right 
DLG subsystem is tuned from V). = 0 to lA = 8 for two 
choices of the right El’s phase: (pr = tt and (pr = 7r/2. 
The corresponding results are displayed in Fig. 4. 

Most notably, at W = 0, a finite dc Josephson cur¬ 
rent only appears if Ap ^ 0. Otherwise Ap = 0 and 
Ap ^ 0 basically cause the same qualitative behavior. 
The amplitude increases with increasing voltage until it 
reaches its maximum at a junction voltage |IT| = 1.0. 
This coincides with the point, W = 2.0, at which the 
El order parameter in the right DLG system attains its 
largest value. If the voltages grows further the ampli¬ 
tude of the exciting Josephson current vanishes rapidly. 
Unfortunately the current becomes extremely small just 
before the El phase completely breaks down. Therefore, 
the Josephson tunnel current is not very suitable for the 
precise determination of the El phase boundary. Figure 4 
corroborates that a finite phase difference yields a larger 
current amplitude. 

Note that only for small gate voltages {W < 0.4) the 
amplitude of the excitonic Josephson tunnel current will 
be of the same order of magnitude as the amplitude of 
a current due to tunneling excitons. For higher voltages 
the normal tunnel current is up to two orders of mag¬ 
nitude larger (see the inset in Fig. 4). That is why we 
propose an experimental setup, where the chemical po¬ 
tentials of the two DLG systems differ only slightly. Nev¬ 
ertheless, we believe that the effect can be observed in a 
wide range of W, simply because the Josephson current 
shows a qualitatively different behavior than the normal 
current. For IF = 0, the only current that may appear 
is the dc Josephson current in the case that the phases 
differ. Any finite gate voltage, on the other hand, in¬ 
troduces an ac Josephson current, while the current of 
unbound electrons is a dc current. 

To conclude, the setup proposed might be used to iden¬ 
tify a condensed exciton phase in DLG, which is sub¬ 
jected to an external electric potential, by analyzing a 



FIG. 3. (Color online) Excitonic Josephson current in case of 
(i) two identical DLG subsystems (green triangles, Vi = 1.0 
and Vr = 1.0), (ii) equal left and right external potentials but 
different phases of the El order parameters (magenta squares, 
Vi = 1.0 and Vr = 1.0), (iii) different external potentials with 
El order parameters having the same phase in the left and 
right DLG systems (blue dashed line, Vi = 0.0 and Vr = 1.0), 
and (iv) different phases and different external potentials in 
both DLG subsystems (red continuous line, Vi = 0.0 and Vr = 
1.0). 
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FIG. 4. (Color online) Amplitude of the excitonic Josephson 
current in dependence on the applied junction voltage in the 
case that the phases of the order parameters are equal {Ap = 
0, blue circles) or differ {Ap = —ttI2, red triangles). The 
dashed and solid lines are interpolations to guide the eye. 
The inset shows the normal tunneling of electrons. The solid 
line considers Akr and, for comparison, the dashed line shows 
the current, when Akr- is set to zero. 


Josephson-type effect in a DLG junction device. Pro¬ 
vided the gate potentials of the DLG systems differ, an 
ac Josephson tunnel is observed irrespective of the phase 
relationship of the excitonic order parameters in both 
subsystems. If both DLG subsystems are exposed to the 
same gate voltages but there is a finite phase difference 
between their exciton order parameters, a dc current ap¬ 
pears. Such a finite phase difference suggests a degen¬ 
eracy of the ground state, i.e., a 1/(1) symmetry. This 
symmetry is closely related to off-diagonal long range or- 
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der and only in this situation the exciting insulator does 
represent a genuine exciton condensate^^. 

We like to emphasize that small leakage currents, 
which may arise in an actual experiment, are linked to 
interlayer hopping of electrons and holes or interlayer ex¬ 
change terms due to the Coulomb interaction. These 
terms pin the phase (j)i to a specific value and there¬ 
fore destroy the C/(l) symmetry^®d9^ Hence exciton 
condensation—in a strict sense—cannot occur and the 
excitonic insulator, if present, features a charge-density- 
wave state. Interlayer hopping and exchange terms such 
as Hu (X moreover, enforce a finite 

Aki at all temperatures and therefore prevent a true 
phase transition^®. 

The present study should be considered as a first step 
towards a theoretical modeling of the excitonic Joseph- 
son effect. The mean-field approach used is certainly 


a crude approximation and any future (more detailed) 
analysis should rely on more elaborated methods that 
take into account fluctuation and correlation effects. Also 
the treatment of the tunnel junction could be improved, 
e.g., by using a more realistic (material specific) tun¬ 
nel matrix element and including pair-breaking effects 
within the barrier. In double bilayer graphene, respec¬ 
tively, double few-layer graphene systems the mobility of 
the electrons and holes is significantly reduced compared 
with double-layer graphene. This reduction strengthens 
electron-hole correlations and corroborates an El forma¬ 
tion at relatively high temperatures®^’®®. For these sys¬ 
tems the theory presented above has to be adapted, but 
the basic idea and the scenario worked out should re¬ 
main valid. Work along these lines will definitely improve 
our understanding of the fascinating exciton condensa¬ 
tion phenomenon. 

This work was supported by the DFG through the spe¬ 
cial research program SFB 652, project B5. 
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